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Abstract 

Using the subdynamical kinetic equation for an open quantum system, a 
formulation is presented for performing decoherence-free (DF) quantum com- 
puting in Rigged Liouville Space (RLS). Three types of interactions were con- 
sidered, and in each case, stationary and evolutionary states were evaluated 
for DF behavior in both the total space and the projected subspace. Projected 
subspaces were found using the subdynamics kinetic equation. It was shown 
that although the total space may be decoherent, the subspace can be DF. 
In the projected subspace, the evolution of the density operator may be time 
asymmetric. Hence, a formulation for performing quantum computing in RLS 
or rigged Hilbert space (RHS) was proposed, and a quantum ControUed-Not 
Logical gate with corresponding operations in RLS (RHS) was constructed. 
A generalized quantum Turing machine in RHS was also discussed. 
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Decoherence, Open System 
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I. INTRODUCTION 



Although an ideal quantum computing system is an isolated quantum system whose evo- 
lution is described by a unitary operator with time reversal symmetry, in practice, it is often 
not practical to isolate a real quantum computing system from its environment. The inter- 
action between a real system and its environment introduces decoherence which destroys the 
superpositions of qubits that enable quantum logical operations to be validated'^'^l. Several 
authors have formulated theories for decoherence-free (DF) subspaces in which quantum 
computing can be performed. Theory is developed starting from a master equation (such as 
the Lindblad equation) for an open system within the powerful semigroup approach. This 
shows that such DF subspaces do indeed exist, allowing logical qubits to be encoded and 
not to decohere^^'^l under the Born-Markov approximation and with restrictions on the type 
of decoherence (i.e., symmetric or collective decoherence). This motivated us to identify 
appropriately designed open quantum computing systems useful for canceling the effects of 
intrinsic decoherence. Such systems would thus behave as ideal quantum computing sys- 
tems in the appropriate subspaces without introducing approximations or restrictions on 
the type of decoherence. In such open systems, self-adjoint operators and unitary evolution 
groups are not intrinsically necessary to govern or operate quantum computation. Quan- 
tum computation can then be performed in a more general functional space which permits 
irreversibility, such as Rigged Hilbert Space (RHS) or Rigged Liouville Space (RLS), rather 
than just Hilbert or Liouville Space. These are kinds of triplet space structure, such as the 
dense subspace and its topological dual space rigged to Hilbert space, $ C ?Y C (RHS) 
or dense space and its topological dual space rigged to Liouville space, ^ <^ ^ G 7i ®TC C. 
x$ (g) (RLS)[^'^l These generalized functional spaces were introduced in physics in or- 
der to provide insight into singular entities like Dirac's delta functions and formulations of 
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quantum mechanicsl^"^! . Since then, they have been intensively developed for describing 
irreversible process, non-integrable systems and chaotic dynamical systems. The latter in- 
cludes RHS of the Hardy class used to formulate the generalized eigenvectors and complex 
spectral decomposition of the LiouviUian for the Priedrichs model'^"^^' and for the general 
spectral decomposition of chaotic maps^^l . Such an approach can also be used to formulate 
open quantum systems in a suitable functional space, revealing non- unitary evolution semi- 
group irreversibility. It may also be suitable for describing quantum logical operations in 
open quantum systems subject to decoherence. 

Based on the above concepts, we propose a type of projected subspace for open quantum 
computing systems, to perform DF quantum logical operations for three different types of 
interactions, and investigate quantum logical operations in such subspaces based on the 
kinetic equation introduced from subdynamics'^'^^l . Compared with recent publications on 
DF subspaces, our approach is general and not restricted to the Born-Markov approximation. 
In section II, a subdynamic formulation and kinetic equation for an open system is derived 
using a simpler approach than previously published references^^'^^l. In sections III, IV and 
V, three types of coherence / DF conditions for the projected subspaces or total spaces 
are analyzed under respective interactions. The role of eigenvalues and eigenvectors of the 
Hamiltonian in the decoherent / DF processes of the systems are discussed. In sections VI 
and VII, a formulation for quantum logical operations in the projected space is investigated, 
and a formulation of quantum computation in RLS (RHS) and a generalized Turing machine 
are proposed. 

II. SUBDYNAMICAL KINETIC EQUATION 

In general, consider an open quantum computing system S, which may consist of real 
spins, pseudospins, or number states of a quantized field (such as photons), coupled with an 
infinite or finite environment R, and where the Hamiltonian operator for the total system is 
given by 
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and Hq [t) = Hs [t) + Hji, Hs {t), and Hi (t) correspond to the Hamiltonian of 5', R 
and the interaction between S and R. A is a couphng constant for this interaction. The 
corresponding Liouville equation for the total density operator p {t) is given by 

ij^p (t) = [Ho (t) + XH, (t) , p (t)] = (Lo (t) + ALi (t)) p (t) . 

Choosing the time-independent eigenprojectors of Lq (t) as and with Qu + Pu = 1, 
respectively, then the eigenprojectors (t) for the total Liouvillian L (t) can be written in 
terms of the Heisenberg equation as i^i^u (t) = [L (t) , H^, (t)]. Hi, (t) satisfy the usual prop- 
erties of projection operators, such as J2u^i> (t) = 1; (0 = 1; (t) (^) = ^u^i, and is 

analytic with respect to A: lim^ >o^u (t) — > P^. From the definition of the eigenprojectors 

11;^ {t) we have QyL (t) {Qy + Py) liy {t) = ZyaQu^u (t) {Zya IS Si correspouding eigenvalue), 
giving 

QyIiy{t)=Cy{t)ny{t), (l) 

where the creation operator is defined by Cyif) =Y^ Cy {Zya,t) PyHyPu with Cy{z,t) = 

a 

'^Zq~JjJ)q~QuL it) Py, for z G complex plane C. In the same way, the destruction operator 
Dy{t) is defined by 

liy (t) Qy = Uy (t) Dy (t) . (2) 

Notice that, since Uy is the projector, therefore Qy ^ Cy{t) ^ Dy (t) in Eqs.([l|) and (||). 
Here Cy {z, t) is the collision operator'^^'^^1 used in non-equilibrium statistical mechanics. 
The operator Cy (t) creates the Qy-pa,Tt of 11;^ (t) from Py and the operator Dy (t) provides 
the opposite operation of the creation operator since Cy (t) = QyCy (t) Py and Dy (t) = 
PyDy {t) Qy. This allows the construction of kinetic equations for each P^-component of 
Uy (t) by 
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p(t) + zp,n,(t)^p(t) (3) 



= [L (t) , (t)] p (t) + p,n, (t) L (t) p (t) 
= p,L (t) (p. + g. (t))n, (t)p(t) 
= p,L (t) (p, + a(t))P.n, (t)p(t) 
= e, (t) (p,n, (t)p(t)), 

where {t) = PyL (t) Pi, + P^L {t) {t) Pi,. A kinetic equation for all of Py is given 

4 ^ (^'^n. (t) p (t)) = E 0. (^) E (^.n, (t) p (t)) (4) 

= e (t) (t) , 

where we define the total intermediate operator as 6 (t) = (^) and the total projec- 

tor as pP™^' (t) = Y.uf^J°^ (t) = Y^yPv^v {t)pit). We refer to Eq.(D as the subdynamical 
kinetic equation, and the operator O {t) as the generator of the subdynamical kinetic equa- 
tion. The second-order approximation for each (t) corresponds to the Bolzmann, Pauli, 
and Fokker-Planck equations of kinetic theory and Brownian motion'^'^^^. The creation and 
destruction operators can be obtained from the basic operator equations within the subdy- 
namics formulation t^'^^l, 

z^a (t) + [Lo, a m = a (a - g.) l, (p, + a) , (5) 

z^D, (t) + [Lo, (t)] = A (P. + D,) Li (g, - D,) . (6) 

The construction of a kinetic equation for the whole system can been accomplished by solving 
the equations for the Ci, (t), or D^, (t), with retarded or advanced integrations; for example, 

a (t) = -ijj^ drfe-*/>°(*')"*'A (a (r) - Q,) U (r) (P, + (r)) fe*/' ^"(*')'^*', 

(7) 

where T is the Dyson time-ordered operator. A similar approach can be used for the de- 
struction operator. Neglecting higher terms {Cl (t)), {D^ (t)), Eqs.(|^) and (||) become 
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i^^D, (t) + [Lo, D, (t)] = -A (P, + D,) L,Q,. 



(8) 
(9) 



The exact solutions for Eqs.(||) and @ are: 



t 

itoo 



Q,L (r) P,e-'^''^°^''"c/r 



JP^LoPut 



(10) 



and 



/±oo 



i r Q^L{T')Q^dT' 



dr 



i Jl^ Q^L(T')Q^dT' 



e J± 



(11) 



The general evolution formulae for the density operator can be deduced from Eq.(0). 
Indeed, since 



and 



p,n, (t) p (t) = (p, + it) a {t))-' (p, + (t)) p (t) 

= (p, + (t) a {t)r' {Pu + 1^. (t)) E (^-' + c.' (t)) p.' (t) 

v' 

= 6,,, (p, + it) a {t))-^ (P.' + (t) a' (0) (0 

V 

= ^p^r (t) , 

where Q is the similarity operator defined as 

V 

Eq.(ll) can be understood as 

= (t) p,n, (t) p (t) 
= e. (t) //;°^- (t) 

= (p.LP, + p.L.a (t) Pu) p^r' {t) , 



(12) 



(13) 



(14) 



(15) 
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which gives the formal solution for fPJ°^ (t) , 

= fe-'^^o "'''"^''^ (P, + (to) a (to))"' {Pu + (to)) (to) . 

Eg. (p!6D , for each P^^- component of Ilu{t) is related to the dynamics of p^^"^ {t) and in- 
cludes the master equation as the second-order approximation of (t) with respect to 
weak coupling'^^1 . The Liouville equation for the reduced projected density of the quantum 
computing subsystem 5" can be obtained by taking the part trace Ttr of Eq. ([l6|) . It can 
be shown that the Liouvillian L (t) and the total intermediate operator O (t) hold a sim- 
ilarity relation (i.e., in the time-dependent case, L{t) = (i^ + n{t)Q{t))n-^{t)). In 
the time independent case, L = QQQ~^, which can replace the Hamiltonian and relevant 
projections'^^'. Therefore, the eigenvectors of in the time-independent case, can be trans- 
formed to L with the same structure of eigenvalues for B and L. Thus Eq.(|15D or (^) is 
the general kinetic equation with significant physical meaning, and can exactly describe the 
dynamics of the open quantum computing system S without invoking any approximations. 
Hence, it provides the starting point for the discussions below. 

III. DECOHERENCE AND DF INTERACTIONS IN TOTAL SPACE FOR A 

DIAGONAL INTERACTION 

Using Eq.(pr5D or (^, it is aparent that the second term on the right tends to zero 
corresponding to a complete DF subspace for the bipartite system S + R. This suggests 
that one should try to identify bipartite systems which, under "diagonal interaction" with 
the environment, still permit the second term to be zero and remain coherent. Here the 
term "diagonal interaction" is taken to mean that the eigenf unctions of the interaction 
part of the Hamiltonian are still the original product states of the unperturbed part of the 
Hamiltonian'^^' . Indeed, if the Hamiltonian of an open system has a diagonal interaction 
form, then the interaction term of Eq.(|r^) or (Q), should equal zero and a DF condition exist 
for the stationary states of the system. 
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To illustrate this point, consider a two-level atom interacting with a single cavity mode. 
Although this model does not describe decoherence from a Bosonic bath (since it only deals 
with one bosonic degree of freedom) it is sufficient to show the underlying principles. In the 
case of an off-resonant interaction, the effective Hamiltonian can be expressed ast^^'^"^! 

H = ujqg^ + ua^ a + ga^ a\2) {2\ , (17) 

where cr^ = |2) (2| — |1) cuq and u are angular frequencies, g is the coupling number, and 
a"*", a are creation and annihilation operators for the cavity mode, respectively. 

Without exercising the diagonal interaction, ga^a |2) (2|, the spectral decomposition of 
the total Hamiltonian for this system can be expressed as 

^0 = E E ((-!)' ^0 + nuj) \j (1) ® n (2)) (2) ® j (1) | (18) 

n j=\ 

with 

2 

1^) = E E (2) ® j (1)1 i^) \ j (1) ® ^ (2)) , for any |?/)) of this system. (19) 

n j=l 

Exercising the diagonal interaction, ga^a\2) (2|, the spectral decomposition of the total 
Hamiltonian is 

^ = E E ((-1)' ^o + nu + ^-^S^ng^ \j (1) ® n (2)) {n (2) ® j (1)| . (20) 

It is therefore apparent that the eigenvectors of this Hamiltonian H are still the origi- 
nal eigenvectors of Hq although the eigenvalues change to (— l)'' a'o + nuj + ^^^^^^ng from 
(— !)•' ujq + nuj. H is diagonal with respect to the basis {\j (1)) ® \n (2)) , (n (2)| (g) {j (1)|}, 
and therefore the off-diagonal matrix elements of the corresponding Liouvillian are equal to 
zero, 

PuLi (t) = 0, for V = (j, f- n, n') , (21) 

Eq.(^ yields 

P.Li (t) (t) = 0, (22) 
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where the eigenprojectors defined in Liouville space are 

= \u) {u\ = \\j (z) ® n (k)) {f (i) ® n' {k)\) (1/ (i) (E) n' (k)) {j (2) ®n{k)\\. 

(23) 

Thus Eq. (|T^ remains invariant with or without the diagonal interaction; that is, there is no 
decoherence introduced by the diagonal interaction for the stationary states. But since the 
eigenvalues have been changed by the interaction, the evolution of the density operators are 
still subject to a type of "decoherence" (unitary error) introduced by this change. That is, 
before the interaction the evolution is 

it)), = E e-<^~'y-"+--y {j (t) ® n (k) 1^ (t)) \j (z) n (k)) , (24) 

jn 

while after interaction, the evolution is 

1^ it))a = E e-K^-^)^'^"^"^-^'^"^)* {j (z) ® n (k) 1^ it)) \j {^) ® n {k)) , (25) 

jn 

and generally we have 

IV^W).^I^W)a- (26) 

In other words, Eqs.(P^ and ( PT| ) reveal that the eigenvectors of the total Hamiltonian H 
with or without the diagonal interaction, have the same tensor products, and are all diagonal 
with respect to the set of these tensor products, showing that the diagonal interaction 
does not introduce any decoherence in the stationary states. Eqs.(p4D and (^) show that 
the eigenvalues of the total Hamiltonian are changed by the interaction. In this case, the 
interaction introduces a phase shift which may introduce a type of "decoherence" (unitary 
error) in the evolution of the states, although the fidelity'^'^^ Tr (^^ p (0) p (t) j , remains 1 for 
both cases (i.e., with or without the interaction). 

IV. DF SUBSPACE FOR TRIANGULAR INTERACTION 

The above discussion for a diagonal interaction, refers to a rather ideal case. Another 
less ideal DF situation is now discussed - it is the "triangular interaction" . This is defined 
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as when the off-diagonal matrix elements (of the interaction part of the total Hamiltonian) 
are up or lower triangular. For example, if the interaction part of the Hamiltonian can be 
expressed in a triangular form using the previous model 



if (1,2) = g (a^aa'' + a a ) , 



(27) 



then the subdynamical kinetic equation (jl^) gives 



i^p. (t) = [PuLsPuP. {t) + PuLbPu + PuLiQ.C, (t) P,] (t) 
= {Ls + Lr) P,p, it) , 

where from the triangular property of the matrix elements we have 

(-1)^ gn, i = j,k = n 

{i®k\H (1,2) \j®n) -- 



(28) 



gVn ~ 1) k = n—l,i = j — l ■ 
0, otherwise 

This causes the second term in Eq.(|T^) to be zero, since the Ci, is a function of QyLiPy from 
subdynamics theoryt^'"^! : 



PuLlQuCuPu — 0. 



(29) 



Eqs.([28|) and ( PU] ) show that the density operator p,y (t) is DF under the triangular interaction 
in the subspaces. That is, with or without loading the triangular interaction, the spectral 
decomposition for the intermediate operator is invariant. 



e = ^P,LoP.= E E {-iyng\jn,fn'){fn',j,n\, 



(30) 



so that an arbitrary expanded projected density operator and its evolution remains invariant 
with or without the interaction, 

2 2 



(31) 



and 



°^°Mt)) = E E k™Mt)) = E E e-(-^)^"^*(/n',j,n|//-Mt))b^,jV)- 



(32) 
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This shows that no decoherence is introduced by the triangular interaction. However, this 
does not mean that the total density operator p (t) is DF in the total space. Indeed, taking 
into account Eq.(|T4|) and using the similarity operator Q on Eg. (pTSP , gives 

i-Qp, (t) = i—p (t) = QQ.Q-'Qp, (t) = QQ,Q-'p (t) , (33) 

and yields 

L = Y, ^Qu^'^ (34) 

V 

= ^ (z/| Lo \u) (P, + a) \u) (z/| (P, + D,) 

= E j: i-iy 9 {n - n') + ((H + EH^.I^')(^'l), 

n/n' j=l \ u'<u J \ u'<u J 

where C^, |z/) can be obtained from the recurrence formula (|^) or ([T0|), 

W\ a (35) 



1 



((-1)^' n - {-lY n') - ((-1)* m - (-1)'' m') 

X [(zm, i'm'\ Li \jn,j'n') + (Im, 2m'| Li |2 (m + 1) , 1 (m' + 1)) {21, 11' \ Cjnj'n' \jn,j'n')] 



with 



{im,i'm'\Li\kl,k'l') (36) 
{im\ Hi \kl) 5k'i'5vm' - {k'l'\ Hi \i'm!) SikSmi 



— gyjm— lSi(^k-l)Sm{l-l)Sk'i'Sl'm' — g\^l' — ^^k' {i' -1)^1' {m' -l)SikSml- 

Hence the eigenvectors of the total Hamiltonian can be formally written as 

M = W)+j:i'^'\CAt)W)W'), (37) 

u'<u 

and 

M = i^\+ E i^\Du{t)W'){iy'\. (38) 

u'<u 

11 



This shows that the eigenvectors of the total Hamiltonian are a combination of the eigen- 
vectors of the intermediate operator 0, and the total system is still subject to decoherence. 
Therefore, it can be concluded that the system, under a triangular interaction in each pro- 
jected Pj, subspace, is DF for both stationary and evolutionary states, but suffers from 
decoherence in the total space. 



V. DF AND DECOHERENT SUBSPACES FOR A GENERAL INTERACTION 

Since the triangular interaction does not introduce decoherence in the corresponding 
subdynamical projected subspaces, this suggests that there should exist a more general 
class of DF subspaces. In fact, as mentioned in the introduction, the existence of such 
DF subspaces has been shown by several approaches, such as projection onto a symmetric 
subspace of multiple copies of a quantum computing system t^l. In this treatment, we focus 
on providing an argument based on a subdynamics approach. Generally, it is assumed that 
the interaction parts of the Hamiltonian for the previous model are non-diagonal or non- 
triangular in form; moreover, for practical use (i.e., to construct a quantum logical gate), a 
model is considered consisting of two atoms with two-levels simultaneously interacting with 
a single optical mode^^^'^^' and influenced by an environment consisting of a set of infinite 
Harmonic oscillators. The Hamiltonian is given by 

H^Ho + XHi (39) 

with 

2 

Ho^Yl [^j^j + 9 + aaf)] + '^a+a + ^ ^kbth, (40) 

j=l k 

Hi = Y.jl9k [hi + h,) (a- + a+) , (41) 

k 3=1 

where 6^, hk are creation and annihilation operators of the oscillator. The eigenval- 
ues of the equation for Hq are given by first writing the matrix expression of Hq in 
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Hilbert space {|± (1) , ± (2) , n, {rife})} and then diagonalizing it. Using the matrix ele- 
ments (± (1) , ± (2) , n', {n'i^}\ Hq |± (1) , ± (2) , n, {n^}), the matrix expression for Hq in the 
subspace {^++, ip-^l^^^^, (^++1,^} , 

(42) 



<1> =| + (l),+(2),n,K}) 



Vnti,k) =|-(l),+(2),n+l,K}), 
^n+i,k) =| + (l),-(2),n+l,K}), 



is given by 



gVn + 1 (n + l)a; + Efe^fe^fc 



^ gVn + 1 (n + l)uj + Efe rikUJk ^ 

Following the eigenvalue equation, -f^o I/) = ^If), gives a characteristic equation, 



0, 



a 


7 


7 


7 







7 
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its solutions give the eigenvalues as: 



(43) 



(44) 



^n,k — 



1 



= 2^ + 2^' + 2 V - 2a6 + a' + 87^), 



2 2 2 
and the corresponding eigenvectors are 



4+1,. = ,:b+-a--^{b^-2ab + a^ + 87^), 



(45) 



/n,fc ) — Vn+l,k) + 



n+l,k 



|6 - |a + i^(62-2a6 + 0^ + 872 
7 



(46) 



n+l,k 



i6 -Ig- 1^(62 _2a 6 + a2 + 872) 



7 



^n+l,kj 
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where we denote a = "^^^"^ +72,0' + J2k''^k^k — b = (72 + 1)0; + J2k''^k<-^k — ^ and 7 



From these eigenvectors, the eigenprojectors Pnm,kj the corresponding unperturbed 
Liouvilhan Lq are defined as 



nm,kj J n,k / \ J m,j J y J m,j / \J n,k 



(47) 



nmjk J yrnm.jk 



with 

^nm,kj nm,kj' 

Without the interaction, the spectral decomposition for Lq is given by 



(48) 



it,=F nm,kj 

= E E (4 

it,=F nm,kj 







^nm,jk^ 


Tnm,jk 



(49) 



k ^m,,j J nm,kj' 



The corresponding spectral decomposition of the intermediate operator 0, in terms of sub- 
dynamics theory, is the same as that for Lq, i.e.. 



^ / y / y nm,kj 0^ nm,kj / , / , \^n,k m,j 
it,=F nm,kj iiT nm,kj 



nmJk J yrnmjk 



(50) 



Now considering interactions occurring to this system, the Liouvillian becomes 

L = Lo + ALi 

with 



(51) 



/ J yrnmjk 
nm,jk^n'm' ,j'k' 



L, 



n'm' ,j'k' 



nm,jk I Vrn'm'jj'k' 



(52) 



Then, in terms of the subdynamical equation (^) the operator B can be expressed as 

/ y / y y nnijkj 0-^ nm,kj ' nm,kj ^^nin,kj nm,jk nm.kj J ' 
it,=F nm,kj 

and the corresponding spectral decomposition is 
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(53) 



e 

±,=F nm,kj 



(54) 



± 



nm,jk 



nm,jk 



nm,jk 



nm,jk 



This again means that the eigenvectors of G, <Pnmjk^ with or without couphng arc the same, 
while the corresponding eigenvalues change to {(ej^ - e^j) + ^iCnmjfc I'^nmjfc)} 

(with the interaction) from (e^^^ — e^j) (without the interaction). Using the set of the 







Tnm,jk J ' 


yrnmjk 



arbitrary state of the system ff^°^^ in the projected subspace space, the stationary projected 
density operator remains invariant with or without the interaction. 



r J / y / J nm,kj nniykjr J / j / j yrnmjk nm,kj nm,kj r 
± nm,kj ± nm,kj 



nm,jk J ' 



(55) 



while the evolution of the projected density operator changes, upon loading the interaction, 
to 



= E E \Pnr^,k,^tZ,kjPit 
± nm,kj 

± nm,kj 



(56) 



k I ^nm,j 



nm,jk 



P±T TT±=F ^ fn 
nm,kj nm,kjr 

yrnmjk nm,kj nm,kjr \ ! J rnmjkj ' 



nmjk J 



± nm,kj 

This is similar to the case for the diagonal interaction, but it takes place in the projected 
subspace, showing that no decoherence occurs to the stationary states, but the evolution- 
ary states are subject to a phase shift introduced by the change in the eigenvalues in the 
subdynamical projected subspace. This phase shift induces a type of "decoherence" (like a 
unitary error) to the evolutionary states. Even so, the fidelity of the mixed states in the 
subspace can be shown to equal 1, 
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F(t) 



(57) 



Tr\l pP''°i (0) re"*/®(*')'^*V^™^ (0) re'/®^*'^''*' 
Tr 



it nm,kj / \ it nm,kj 



Tr 



\ 



it nm,kj 



nm,kj 



^ l\ n,k m,,j J ' \^ nm,jk \ nm,jk \ ^ nm,jk J J 



y nm,kj^ 



Tr 



\ 



lnm,kj 



nm,kj 



it nm,kj 

Tr U pP-°i {t^) pP-°i {to) 



where {f^^"-')^ kj matrix element of the density operator pP^'°^ (to) with respect to 



6^^ ) 




rnm,jk J 1 


yrnmjk 



Since the eigenprojectors remain invariant, the initial density operator remains the same 
with or without the interaction under the assumption of decoupling of initial states. Eq.(p7D 
shows that the constructed projected subspace is DF! Whatever the subspace, the total 
density operator p is decoherent in the total space. This suggests a universal procedure for 
constructing a DF subdynamical projected subspace by choosing a set of one-dimensional 
eigenprojectors of the unperturbed Liouville operator Lq (t) through the kinetic equation 
(|15|) or (H). In this constructed subspace, the governing equation is the subdynamical kinetic 



equation and the interaction part of the total Liouvillian cannot change the diagonal property 
of the intermediate operator, if each of Pi, in the complete set is chosen as a one-dimensional 
projector by construction (note if Pj, is over one-dimensional the intermediate operator 
generally is not diagonal). 

The above discussions are summarized below in tabular form: 
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Interaction Type 


Total Space 


Projected Subspace 




St 


Ft 


St 


Ft 


Diagonal Interaction 


DF 


DF 


DF 


PE 


Triangular Interaction 


DF 


DF 


DF 


DF 


Diagonal 


D 


D 


DF 


PE 



where St and Et refer to stationary states and evolutionary states, and PE, D and DF 
refer to a phase error, decoherent and decoherence-free, respectively. The table shows that 
it is possible to construct a DF subdynamical projected subspace for an open quantum com- 
puting system with three types of interaction occurring with the environment. This allows 
quantum logical computations to be performed by using the projected density operator on 
the corresponding subspaces, which may be beyond Hilbert space. For example, these can 
include rigged Hilbert spaces. To clarify this point, a model of two two-level atoms is dis- 
cussed in detail in the following subsection. A second-order calculation is presented, but 
the approximation does not restrict the validity of the result, since there is no problem in 
allowing high order calculations to be performed based on the same methodology. 



VI. QUANTUM LOGICAL OPERATION IN SUBSPACES 



The expected evolution of the projected density operator, in the projected subspaces, 
can be formally solved from Eq.(^), 



it^ nmjk 



(58) 



EE' 

di=F nmjk 



nm,kj ^ 



yrnmjk rnm,kj \ / J rnm,jk J 



where the difference in energy E^^^Jf is defined as E^^j^j + AEnm,kj = ~ ^mj + 
A (^<Pnmjk ^i^nmjk ^nmjk^^ ^he projected density operator is determined from p(t). 



p±±TTY\±±TT r,(f) = p-^^TT (f) 
nm,kj nm,kj r y'^J rnm,kj \ J 



(59) 



y nm,kj ~ nm,kj ^nm,kj J y nm.kj ' nm,kj J I \^ J 
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Here D^^Jf is given by Eq.(|TT]), 



nm,kj 

— Q AP nm.jk nm.jk 



S/ p ^ nm.jk ^ nm,j k 



(60) 



_ \ p±±TT r /n±±TT i 

nm,jk ^l^nrnjk ^±±=p=p ^ ^±±=p=ji 



Taking into account the Born approximation, the second order term D^^j^fj^^ is 



D 



[2]±±TT 
nm,jk 



-A 



j2 

0: 



±±=FT 
nrnjk 



A 


Tnm,jk 




'rn"m",j"k" ) 








- 


\^ n.k m,j 


) ±zO 



,jk ) yrn"m",j"k" 



m,jk 



^n"m",j"k" 



+ A^ , 



(61) 



which gives the second order projected density operator as 



rnm,kj \ J nm,kj ' ni 



im,kj 



Pit) 



(62) 



where | L? = g, {^U + l)j + ^ [j - I) j) - Qk {^{k + l)k + ^{k-l)k). 

On the other hand, E^^JJ can be approximately found by calculating C^^jk from 
Eq.(p!oD. For example, the first order cl^^^^j^^ is 



nm,jk 



-A- 



1 



^nm,jk ^nm.jk \^n,k ^m,j J 
1 



f)±±TTT p±±T=F 
^nm,jk -'^1-' nm,jk 



(63) 



±T n"m"j"k" {^n",k" ^m'',j" 



+ ± zO 



^ (^n'W'j-"fc" I -^1 0nmj?) ^n'W^J''fc'') (^nmj? + (A) , 

which then gives the second order term E^^j^^^j^^ , 



E. 



[2]±±TT 
nm,kj 



S-^;S)-A'E E 





r^nm,jk 




r'nmjk 1 


\^n",k" 


^-m",j" J 


- 


^ n.k ^m,j J 





(64) 



y'n,k m,j J 



AE E 

it^ n'm',j'k' 



g, (^(j + 1) J + y (j - 1) j) - g, {^{k + l)k + ^{k- 1) fc) 
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and hence giving the second order evolution of p?"'°-?Pl (t), 



p.om (,)) ^Y.Y. kitr^ (^)) = E E e-^""-^""* kSr^ (0)) \&) ■ 

±=F nmjk ±=F nmjk 

(65) 

Considering this evolution, a quantum XOR operator for the system may be constructed by 
applying a sequence of operations'^^'^^l which is related the swap operator, 



(66) 



Suppose there is no coupling interaction, the ideal swap operator is then given by adjusting 
the coupling time tgw 



±^ nmjk 
±T nmjk 



(67) 



4^nm,jk ) 









and loading the interaction, the nonideal action of the swap operator is given by the evolution 
operator (since the energy shift induced an error), i.e.. 



±=F nmjk 



nmjk J \rnm,jk 



To cancel the decoherence, we have to have 



EEe- 

±T nmjk 







T^nmjk J 


^Tnmjk 



(68) 



(69) 



which gives 



sw — EO . ; 



(70) 



where we have assumed that the distribution of energy is homogeneous, ^ = constant, to 
determine the universal time At. For example, if the second order nonideal swap operator 
is 
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it^ nmjk 



nmJk J yrnmjk 

then the second order At[^], is 



AtE (72) 

3j (v^(i + lji+V^(Fl)i) -9k (^^y{k+l)k+^y{k-l)k^ 



tsw^'l2±^'l2n'm' ,j'k' ^ ^ ±± ip^ \ / ±± 



This shows that: (1) the change of the eigenvalues after loading the couphng interaction, 
introduces a phase error in the swap operator, (2) the error can easily be cancelled by 
adjusting the coupling time (At^^) based on tsw, since the eigenvectors are invariant. This 
demonstrates an important advantage of performing quantum Logical operations in the 
projected subspace. 

Moreover, from Eqs. (p!oD , (PD and (|6^) it is evident that eigenvalue -E'^^JJ be com- 



plex (in fact, this model is a type of Friedrichs model, possessing a complex spectrum in 
the sense of extending Hilbert space'^^"^''!). This means that the time evolution of Pnm%J 
in Eq.(|58D may be asymmetric. Thus, a Rigged Liouville Space (RLS) formulation can be 
chosen'^'^1 to describe the subsystem, since a Hilbert space formulation provides an inade- 
quate description. The evolution of the projected density operator and its adjoint operator 
are expressed as 

E E (^)) = E E e-^^-""* \pi^%j (0)) l^^^S^) , (73) 

it=F nmjk it=F nmjk 

E E {P^mjj it)\ = E E {p^mJJ (o)| , 

it^ nmjk it^ nmjk 

where Pnm%J i^) exist in the test space ^sr ® ^sr, which is a dense subspace of the 
Liouville Space, representing the physical states which can be prepared in an actual exper- 
iment. Its adjoint p^j^ ^J (t) lies in the dual space ^^sr ® ^sr^ representing a procedure 
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that associates with each state a number, while preserving the hnear structure which results 
from the superposition principle. This is a RLS structure which facilitates describing irre- 
versible processes like decoherence and dissipation due to interaction with the environment. 
This poses an interesting question as to what is new in this formulation for an open quan- 
tum system, if it performs quantum logical operations in a RLS (RHS). For convenience 



we assume that 



Pnm^kj if))i {Pnmjtj (^) biorthouormalized with respect to each other, 
and the evolution of Pnm%J i^) i^^ Eq. (|58|) ) in the projected subspace permits the system to 
perform quantum computing in a RLS (RHS). 

Indeed, if we construct a general quantum logical operator Q given by some function / 
of the biorthonormal states | Pnm,kj (^)) ; {Pnm,kj (t) | in a RLS as 

Qg {\PnZ,kj (t))) = 9 ^PnZk, W)) , (74) 



and 



(75) 



(76) 



2 — / ( Pnm,kj (O) {Pnm,kj (t) 

Then the states or observables are g [ Pnm,kj i^))) ^ g' { p'nm,kj (t))) e $5 ® $s, 
h (^{p^j^kj (t) h' [(^Pnm,kj (^) ) ^ ^^s®^^-, and the space is closed to the dual pair spaces 
$5 ® $5 and $5 (g) $5 . The RLS structure is $5 $5 C 7^ C $5 $5 . For example, 
a quantum universal gate, such as a two qubit based quantum Controlled-Not logical gate 
(XOR gate) in RLS, can generally be constructed from a combination of the projection 
operators as: 

CM = |00) (00| + |01) (0T| + 1 10) (lT| + 111) (10|, (77) 
which induces the corresponding quantum logical operations, 

CX |00) = |00) , CAf |01) = |01) , (78) 
CAf\10) = |ll),CAr|ll) = |10), 
(OOICAT = (00 1, {01\CAf = (01 1, 
(J6\CAf = (Til, (Jl\CAf = (To I, 
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and demonstrates that the quantum logical operations are closed in RHS. The property of 
closure for quantum logical operators in RLS (RHS) is important for practical applications. 
This shows that quantum logical operations in RLS (RHS) inherit the essential properties of 
those in Liouville (Hilbert) space, if the biorthornomal property of generalized eigenstates 
still hold. 

One may still argue that the advantages of the above approach for quantum logical opera- 
tions in RLS (RHS), for an open computing system S, are not evident. We would argue that 
the Liouville (Hilbert) space approach generally fails to describe open quantum computing 
systems which undergo irreversible processes. These include semigroup evolution and non- 
self-adjoint Hamiltonians of subsystems, far from equilibrium. This is something that the 
RLS (RHS) formulation is well suited to handhng. To construct quantum logical operations 
in projected subspaces, the influence of decoherence induced by changes in the eigenvalues, 
may easily be cancelled since the eigenvectors remain invariant. Moreover, in this space 
the evolution of the states are permitted to be time asymmetric, providing a framework for 
describing the irreversibility of practical open systems. This irreversibility does not change 
quantum reversible logical operations to quantum irreversible logical operation in the quan- 
tum universal Controlled- Not logical gate. To appreciate this, one must distinguish between 
irreversibility of a quantum logical operation, introduced by the structure of logical gate 
and irreversibility of the process, induced by interactions with the environment. Reversible 
computation means reversible logical operations on the structure of the logical gate. In this 
sense, quantum computing in RLS (RHS) does not adversely impact on reversible quantum 
logical operations and permits computing any reversible function, although irreversible pro- 
cesses exist. This is further clarified below in relation to construction of a quantum Turing 
machine in RHS. 
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VII. CONSTRUCTION OF A TURING MACHINE IN RHS 



Consider, in general, an open system with N components (atoms) each of which can 
have two states. We design a quantum Turing machine consisting of N = n + 1 pseudospins 
or levels |±(j)), j — h,l,2,--- ,n, in a 2"+^ dimensional test space $ and ^± (j) , j — 
h,l,2, - ■ ■ ,n, in the dual space of the test space with biorthnomality and isometry. 



±(J) 



(79) 



± m 



(80) 



where the Turing head is denoted as h, and the Turing tape pseudospins are denoted as 
1, • • • ,n, and ± represent two states or two observables. The quantum Turing machine can 
be considered as an open system in which the Turing tape is a special finite bath coupled to 
the Turing head. The state of the quantum Turing machine lies in a 2""*"^ dimensional 
dense subspace $ of the Hilbert Space and is spanned by {|± (j))}, and its observables (V^ 



arc spanned by <{ ( ± (j) 
given by 



in the dual space The evolution of a pure state l-i/^o) (V'o 



IS 



=c/(j)IV'o)(^o u{j) 



(81) 



where the non-unitary evolution operator U can be expanded by the non-self-adjoint gener- 
ators defined by 



KU) =^oi(j) + Ao(i), 

A, (j) = ^Pol (j) - tPio (j) 
K U) = Pn U) - Poo U) , 



(82) 



with 



Pik{j)^\i{j))Ck{j) 



(83) 
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Following ref . [25] , we restrict ourselves to the Bloch vector A of the Turing head h as 



Xl = Uh A, (/i) ® 1 (1) ® ■ ■ • ® 1 (n) liJh) . 



(84) 



Hence, for 2" biorthonormal initial tape states |(^o) = |='=(1)) ® \^{'^)) ® ••■ ® 



00 =(±H 



±(2] 



®(±(1) 



and arbitrary Turing head h states \Lph) , {fh\, the 



states of Turing machine \iph), (iph can be written as the tensor product states 



(85) 



and hence the Bloch vector of the Turing head (Q) is described by the above tensor product 
states as 



Af (00, 0o) = {Vh ® </)o| Ai (/i) ® 1 (1) ® ■ ■ ■ ® 1 (n) \iph (S) (po) 



(86) 



performing a pure state trajectory on the Bloch circle, 



Aj (0o,0o)' + A^ (00, 0o)' 



1. (87) 

Any initial state of the Turing machine iV'o), ("00 with the Turing head h in the pure state 
\{po (h)), {tpo {h)\ can thus be written 



(88) 



2" 



i=i 

This yields a motion of the Bloch vector of Turing head h induced by the change of the tape 
states, given by 



2" 



(^o,^o) =E«A A^ (00, 0^ 



(89) 
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where aj, bj are expansion coefficients. This shows that the motion of the Bloch vector of 
Turing head (^i/jq, V'o) constructed in RHS can be expressed as 2" — 1 entangled Bloch vec- 
tors of the Turing head (0O)0o) induced by a pure Turing head state \ipo (h)), {(po {h)\, 
and 2" initial tape states |0o), (^o ■ This means that the motion of the Turing head h, 
X'l ^Vo) V'o); can perform many parallel primitive trajectories of the Bloch vector of the Tur- 
ing head h, A^ (^(pi, (plj, identical to the situation of the reversible Turing head constructed 
in the Hilbert Space formulation I^^l . This is because the generahzed basis in RHS remains 
biorthonormal and isometric, corresponding to the orthonormal property, and invariance 
of the inner product for the basis in Hilbert Space. This enables most properties of the 
quantum Turing machine in Hilbert Space to be inherited in RHS, with the exception of the 
reversible property of an ideal process. 



VIII. CONCLUSIONS 

In conclusion, a subdynamic formulation for an open quantum system is presented using 
a simpler approach than discussed previously. Based on the subdynamical kinetic equation 
for an open quantum system, the properties of coherence and decoherence-free conditions 
for three types of interactions were analyzed. A condition for a quantum computing system 
to remain in coherent states in the projected subspaces is found. That is, the interactions 
between the system and its environment should be "diagonal" or "triangular". Moreover, 
one can find universal DF projected subspaces by using the subdynamical kinetic equation 
without restrictions on the type of decoherence and without introducing any approximations, 
although the total space may be decoherent. An implied universal concept in our analysis is 
that the invariant eigenvectors remain a DF structure for any stationary states of the system 
and that changes of the eigenvalues only induce a phase shift. This type of phase shift which 
may introduce a type of phase error for the evolved states does not change the DF property 
of the subspace. The possible phase error can be cancelled, due to the invariance of the 
eigenvectors in the projected subspace under the assumption of a homogeneous distribution 
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of the energy shifts. In this projected subspace, one can construct a quantum XOR operator 
and show that the evolution of the SWAP operator is time asymmetric. Hence, we propose 
a formulation for performing quantum computing in RLS (RHS) by constructing a general 
quantum Controlled- Not logical gate, with corresponding operations in RLS (RHS), and a 
generalized quantum Turing machine in RHS. 
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